This Page Is Inserted by IFW Operations 
and is not a part of the Official Record 

BEST AVAILABLE IMAGES 

Defective images within this document are accurate representations of 
the original documents submitted by the applicant. 

Defects in the images may include (but are not limited to): 

• BLACK BORDERS 

• TEXT CUT OFF AT TOP, BOTTOM OR SIDES 

• FADED TEXT 

• ILLEGIBLE TEXT 

• SKEWED/SLANTED IMAGES 

• COLORED PHOTOS 

. BLACK OR VERY BLACK AND WHITE DARK PHOTOS 

• GRAY SCALE DOCUMENTS 

IMAGES ARE BEST AVAILABLE COPY. 



As rescanning documents will not correct images, 
please do not report the images to the 
Image Problem Mailbox. 



CM 
O 
CM 



(19) 



J 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(12) 



(id EP1 202 119 A1 

EUROPEAN PATENT APPLICATION 



(43) Date of publication: 

02.05.2002 Bulletin 2002/18 



(51) intci ? G03F7/20 



(21) Application number: 01307653.4 

(22) Date of filing: 10.09.2001 



(84) Designated Contracting States: 

AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 
MC NL PT SE TR 
Designated Extension States: 
AL LT LV MK RO SI 

(30) Priority: 12.09.2000 US 231953 P 

(71) Applicant: ASML Masktools Netherlands B.V. 
5503 LA Veldhoven (NL) 



(72) Inventor: Llebchen, Armin 

Hayward, California 94542 (US) 

(74) Representative: Leeming, John Gerard 
J.A. Kemp & Co.. 
14 South Square, 
Gray's Inn 

London WC1 R 5LX (GB) 



(54) 



Method and apparatus for fast aerial image simulation 



(57) The present invention provides a method and 
apparatus for simulating an aerial image projected from 
an optical system, wherein the optical system includes 
a pupil and a mask. In general, the method comprises 
the steps of obtaining parameters for the optical system, 
calculating a kernel based on an orthogonal pupil pro- 
jection of the parameters of the optical system onto a 
basis set ; obtaining parameters of the mask, calculating 
a vector based on an orthogonal mask projection of the 
parameters of the mask onto a basis set, calculating a 
field intensity distribution using the kernel and the vec- 
tor, and obtaining aerial image data from the field inten- 
sity distribution. 



Fig.5. 

Ilium shape 
sigma 



specify discrete grid 
for pupil/illuminator 
sampling 



Discretize illuminator profile 
over tinite grid and sample 
weights w_s 

7 



tabulate array ot 
j[a](k_i) over 1-D pupil 
grid k_ i. 
EQ.Q9), (20) 




Calculate p(a](s J)over 
pupil grid using a zero 
model EQ. (95) 




opd 



Calculate pf>](s_i)over 
pupil gnd usmg diffusion 
model EQ. (t13) 



Calculate p[a](s_i) 
over puDil grid using 
opd model EQ.(85) 



Combine w_s and 
p|a](s_i) EQ. (96) 




< 



CL 

LU 



Pnntcd by Jouvc. 75001 PARIS (PR) 



EP 1 202 119 A1 



Description 



BACKGROUND OF THE INVENTION 



w 



15 



25 



30 



35 



40 



45 



50 



55 



[0001] In lithography, an exposure energy., such as ultraviolet light that is generated in an optical system, is passed 
from an aperture of the system through a mask (or reticle) and onto a target such as a silicon substrate. The mask 
typically may contain opaque and transparent regions formed in a predetermined partem. The exposure energy exposes 
the mask partem, thereby forming an aerial image of the mask. The aerial image is then used to form an image onto 
a layer of resist formed on the target. The resist is then developed for removing either the exposed portions of resist 
for a positive resist or the unexposed portions of resist for a negative resist. This forms a patterned substrate. A mask 
typically may comprise a transparent plate such as fused silica having opaque (chrome) elements on the plate used 
to define a pattern. A radiation source illuminates the mask according to well-known methods. The radiation transmitted 
through the mask and exposure tool projection optics forms a diffraction-limited latent image of the mask features on 
the photoresist. The patterned substrate can then be used in subsequent fabrication processes. In semiconductor 
manufacturing, such a patterned substrate can be used in deposition, etching, or ion implantation processes to form 
integrated circuits having very small features. 

[0002] In a manufacturing process using a lithographic projection apparatus, a pattern in a mask is imaged onto a 
substrate, which is at least partially covered by a layer of radiation-sensitive material (resist). Generally, lithographic 
patterning processes are understood by those who practice the profession. Information regarding exemplary processes 
may be obtained, for example, from the book "Microchip Fabrication: A Practical Guide to Semiconductor Processing- 
Third Edition, by Peter van Zant, McGraw Hill Publishing Co., 1997 ISBN 0-07-067250-4. 

[0003] As the size of lithographically fabricated structures decreases, and the density of the structures increases 
the cost and complexity of designing masks additionally increases. One method of reducing costs of lithographic fab- 
rication is by optimizing the lithographic design with a lithographic simulation stop prior to the actual manufacturing 
step. One specific method of lithographic simulation is drawn to simulating the aerial image of the mask. The aerial 
image is defined as an intensity distribution of light just prior to reaching the resist on a surface of a substrate, when 
the substrate is exposed via the mask in an exposure apparatus. In order to simulate the aerial image, a layout of a 
mask and exposure conditions (non-limiting examples include NA: Numerical Aperture, a (sigma): Partial Coherence 
Factor) of the lithographic apparatus are typically required as input parameters. 

[0004] Lithographic apparatus may employ various types of projection radiation, non-limiting examples of which in- 
clude light, ultra-violet ("UV") radiation (including extreme UV ("EUV"), deep UV ("DUV"), and vacuum UV ("VUV")). 
X-rays, ion beams or electron beams. The following have been considered exemplary exposure sources. Light may 
generally refer to certain mercury emissions, i.e., wavelengths of 550 nm for the Mine, 436 nm for the g-line and 405 
nm for the h-hne. Near-UV or UV generally typically referto other mercury emissions, i.e., 365 nm for the Mine DUV 
generally refers to excimer laser emissions, such as KrF (248nm) and ArF (193nm). VUV may referto excimer laser 
F 2 , i.e., 157nm : Ar 2 , i.e., 126nm, etc. EUV may referto 10-15 nm. This last portion of the electromagnetic spectrum is 
very close to "soft X-rays" but has been named as "EUV", possibly to avoid the bad reputation of X-ray patterning Soft 
X-rays may refer to 1 - 15 nm, which may typically be used in X-ray lithography. 

[0005] Depending on the type of radiation used and the particular design requirements of the apparatus, the projection 
system may be for example, refractive, reflective or catadioptric, and may comprise vitreous components, grazing- 
incidence mirrors, selective multi-layer coatings, magnetic and/or electrostatic field lenses, etc; for simplicity such 
components may be loosely referred to in this text, either singly or collectively, as a "lens". 

[0006] When the resist is exposed by the aerial image, there is an additional variable in that some of the exposure 
light is reflected back by the surface of the substrate, and then absorbed by the resist. Accordingly, not only the resist 
charactenst.es (regarding, for example, refractive index: Dill's A,B,C) but also parameters regarding the characteristics 
of the substrate (e.g., refractive index) should be included as input parameters for simulating the latent image 
[0007] The so-called Hopkins model treats the electric field forming the image typically as a scalar and assumes the 
object being imaged is sufficiently thin so that its effect on the incident field is represented by a multiplicative function 
It is advantageous to perform the image formation analysis in the Fourier domain (frequency space) in order to deal 
with the pupil function of the imaging system rather than the amplitude response function, and with the angular distri- 
bution or "effective source" rather than with the mutual intensity. 

[0008] There arc several computer programs commercially available that calculate aerial images based on the Hop- 
kins model. For example, the University of California at Berkeley, Department of Electrical Engineering and Computer 
Science, Berkeley, Calif., 94720, offers a program called SPLAT. 

[0009] The Hopkins model is used to model the imaging of drawn design features under partially coherent illumination 
A major problem in modeling aerial images under partially coherent illumination is the necessity to superimpose and 
add the effect of each individual illumination source that makes up the partially coherent source. In the Hopkins model 
a two-dimensional by two-dimensional transmission cross-coefficient function ("TCC") is pre-calculated. which captures 
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all the effects of the lithographic projection apparatus, including NA. sigma. etc. As taught, for example in Born & Wolf , 
p. 603. once a TCC is known, systems with partially coherent illumination can be modeled by integrating the TCC over 
the Fourier transform of the transmission function for the geometrical layout feature under illumination. 
[0010] Fundamentally, the TCC is a two-dimensional by two-dimensional correlation function with a continuous set 

5 of arguments. In practice, an assumption can be made that the feature patterns to be imaged are periodic in space. 
For such periodic patterns, the TCC has a large, but discrete, set of arguments. The TCC can then be represented as 
a matrix, with discrete columns and rows. For typical features interesting to the lithographer the size of this matrix is 
tremendous and restricts the scale and size of features that can be simulated. It is a purpose of any simulation algorithm 
based on the Hopkins model to reduce the size of this matrix by approximation, while retaining a reasonable degree 

10 of accuracy. 

[001 1 ] An exemplary projection lithography system is illustrated in FIGS. 1 A and 1 B. In FIG. 1 A, light from illumination 
source 102 is focussed by condenser 104. The condensed light passes through the mask 106, then through the pupil 
108 and onto the substrate 110. As shown in FIG. 1B. substrate 110 may comprise a top a nti- reflective coating 112, 
a resist 118, a bottom anti-reflective coating 117, a top substrate layer 114 and a plurality of other substrate layers. 
15 As indicated in FIG. 1 B, the focal plane may lie within the resist 118. 

[0012] In the past, several numerical techniques have been applied to reduce the size of the TCC to reasonable 
scales. In one instance singular value decomposition has been applied to decompose the TCC into its eigenspectrum, 
sort the resulting eigenvectors in decreasing magnitude of their eigenvalues, and only retain a finite number of eigen- 
vectors in order to approximate the TCC. An exemplary method for optical simulation for the system of FIG. 1A. thai 
uses the Hopkins model, is illustrated in FIG. 2. As illustrated in FIG. 2, at step S202, the lithographic projection ap- 
paratus and mask parameters are input in the system. At step S204, the TCC is approximated by an eigenvector 
diagonalization. At step S206, a two-dimensional table of convolution integrals of the eigenf unctions with a discrete 
set of normalized rectangles is calculated and tabulated. The article, "Fast, Low-Complexity Mask Design", by N. Cobb 
ct aL SPIE Vol. 2440, pgs. 31 3-326, the entire disclosure of which is incorporated by reference, teaches an exemplary 
method for completing step S206. Since the resulting table is on a discrete grid of possible layout-features, to attain 
an acceptable resolution, the table must be large. At step S208, since the TCC has been approximated, the illumination 
system can be modeled, i.e. the aerial image can be simulated by combining each two-dimensional convolution integral 
corresponding to each respective rectangle in a given proximity window. 

[0013] This method has two drawbacks. First, an eigenvector diagonalization is a numerically expensive operation. 
Thus, for a user to change lithographic projection apparatus conditions, such as NA, sigma, illumination type or lens- 
aberration, an expensive recalculation of the TCC approximation that limits the use of the simulation tool is required. 
Second, the eigenspectrum of the TCC is a two-dimensional function. Representing the associated field vectors of the 
geometrical pattern under illumination requires a two-dimensional lookup table that limits the speed of the field calcu- 
lation in the Hopkins algorithm. Specifically, since the field vectors are a two-dimensional table, two pointers are needed 
to access any respective address of data, which increases access time. Further, since the field vectors are a two- 
dimensional table, they are stored in a portion of the cache. Using a look-up table from a DRAM in a present day CPU 
is more expensive than performing multiplication because the cache miss penalty is larger than the time to perform 
the multiplication. 

[0014] Accordingly, there is needed a method and apparatus for simulating a projection lithography system using 
the Hopkins model, that does not use a numerically expensive eigenvector diagonalization operation. What is further 
needed is a method and apparatus for simulating a projection lithography system using the Hopkins model that permits 
a user to change lithographic projection apparatus conditions, such as NA, sigma, illumination type or lens-aberration 
without requiring an expensive recalculation of the TCC approximation. 

45 BRIEF DESCRIPTION OF THE INVENTION 

[001 5] 1 1 is an object of the present invention to provide a method and apparatus for simulating a projection lithography 
system using the Hopkins model that does not use a numerically expensive eigenvector diagonalization operation. 
[0016] It is another object of the present invention to provide a method and apparatus for simulating a projection 
lithography system using the Hopkins model that permits a user to change lithographic projection apparatus conditions 
without requiring an expensive recalculation of the TCC approximation. 

[001 7] It is still another object of the present invention to provide a method and apparatus for predicting the intensity 
field distribution ("aerial image") at the surface and throughout the resist film of a substrate being irradiated in a litho- 
graphic imaging process. 

[0018] In accordance with the foregoing objectives, the present invention writes the TCC as T(q*q). where q and q' 
each are continuous frequencies in two-dimensional Fourier space, q=(q x , q y ), q'=(q' x , q' y ). The present invention ap- 
proximates the TCC extremely well as a bilinear form of a basis function with kernel A[ij]. Furthermore, the present 
invention utilizes a set of orthogonal polynomials. The kernel A[ij] represents a small matrix that efficiently approximates 
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the TCC for a broad range of illumination conditions used in modem lithographic processing. The present invention 
uses a lower number of arithmetic operations as a result of the use of a one-dimensional look-up table as opcosed to 
the two-dimensional look-up table in the prior art. Furthermore, the one-dimensional look-up table achieves a higher 
cache hit ratio than a two-dimensional look-up table of the same resolution, thereby providing a more efficient system. 
[0019] in general, in one embodiment, the invention features a method of simulating an aerial image projected from 
an optical system, the optical system including a pupil and a mask plane, the method comprising the steps of providing 
a mask to the mask plane, obtaining parameters for the optical system, calculating a kernel based on an orthogonal 
pupil projection of the parameters of the optical system onto a basis set, obtaining parameters of the mask, calculating 
a vector based on an orthogonal mask projection of the parameters of the mask onto a basis set. calculating a field 
intensity distribution using the kernel and the vector, and obtaining an aerial image from the field intensity distribution. 
[0020] In one embodiment of the invention, the parameters for the optical system include aberrations. 
[0021] In another embodiment of the invention, the step of calculating a kernel corresponding to the parameters of 
the optical system includes the step of generating a table of indefinite integrals based on a recurrence over a seed 
array of incomplete gamma functions. 

[0022] In yet another embodiment of the invention, the step of calculating a kernel corresponding to the parameters 
ol the optical system further includes the step of tabulating an array of orthogonal pupil projection coefficients corre- 
sponding to respective points in the pupil of the optical system, wherein the optical system is in-focus. More particularly, 
it further comprises the step of combining sample weights of an illuminator profile of the optical system with the array 
ol orthogonal pupil projection coefficients. 

[0023] In still another embodiment of the invention, the step of calculating a kernel corresponding to the parameters 
of the optical system further includes the step of tabulating an array of orthogonal pupil projection coefficients corre- 
sponding to respective points in the pupil of the optical system, wherein the optical system is either not in-focus or has 
aberrations. More particularly, it further comprises the step of combining sample weights of an illuminator profile of the 
optical system with the array of orthogonal pupil projection coefficients. 

[0024] In still yet another embodiment of the invention, the step of calculating a kernel corresponding to the param- 
eters of the optical system further includes the step of: tabulating an array of orthogonal pupil projection coefficients 
corresponding to respective points in the pupil of the optical system, wherein the optical system accounts for effects 
of diffusion of photoactive compounds in the resist. More particularly, it further comprises the step of combining sample 
weights of an illuminator profile of the optical system with the array of orthogonal pupil projection coefficients. 
[0025] In a further embodiment of the invention, the step of calculating a vector corresponding to the parameters of 
the mask further includes the step of specifying a proximity window within the mask for geometric sampling. 
[0026] In still a further embodiment of the invention, the step of calculating a vector corresponding to the parameters 
of the mask further includes the step of decomposing a geometric pattern of the mask into a disjoint set of rectangles 
and tabulating an array of projections of the rectangles within a proximity window. More particularly, the step of calcu- 
lating a vector corresponding to the parameters of the mask further includes the step of correcting the array of projec- 
tions of the rectangles based on the type of mask. 

[0027] In general, in another aspect, the invention features a method for simulating an aerial image projected from 
an optical system, the optical system including a pupil and a mask plane, the method comprising the steps of: providing 
a mask to the mask plane, obtaining parameters for the optical system, obtaining parameters of the mask, projecting 
the components of electric field vectors to an orthogonal polynomial basis as described in further detail below, the 
polynomials including polynomials associated with the parameters of the optical system, and approximating the trans- 
mission cross-correlation function associated with the optical system based on the orthogonal projection of polynomials. 
[0028] In general, in still another aspect, the invention features a simulation device operable to simulate an aerial 
image projected from an optical system, the optical system including a pupil and a mask plane, the simulation device 
comprising a first parameter obtainer for obtaining parameters of the optical system, a first calculator for calculating a 
kernel based on an orthogonal pupil projection of the parameters of the optical system onto a basis set, a second 
parameter obtainer for obtaining parameters of a mask provided to the mask plane, a second calculator for calculating 
a vector based on an orthogonal mask projection of the parameters of the mask onto a basis set, a third calculator for 
calculating a field intensity distribution using the kernel and the vector, and an aerial image obtainer for obtaining aerial 
so image data from the field intensity distribution. 

[0029] In one embodiment of the invention, the first, second, and third calculators are the same calculator. 
[0030] In another embodiment of the invention, the parameters for the optical system include aberrations. 
[0031] In yet another embodiment of the invention, the first calculator is operable to generate a table of indefinite 
integrals based on a fast recurrence over a seed array of incomplete gamma functions. 
55 [0032] In still another embodiment of the invention, the first calculator is operable to tabulate an array of orthogonal 
pupil projection coefficients corresponding to respective points in the pupil of the optical system, wherein the optical 
system is in-focus. 

[0033] In still yet another embodiment of the invention, the first calculator is operable to tabulate an array of orthogonal 
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accompanying drawings. 
BRIEF DESCRIPTION OF THE DRAWINGS 



In the drawings: 



r,^ S o 1 B illUStra,e a conv *ntional optical lithographic system 

» with J, exempt^ imaae P«*«" <™ « optica, system in accordance 

K FIG.' 6 ESS a a n" 1ZS Sod" !o°: oTT 9 3 ^ * P — '-ntion. 

[0046] FIG. 7 illustrates .TJZ^^^^S^ ! T*!* aCCOrdanCe with the P re ^t invention. 

S3 FIG. 9 X^ZZZZSZr i^ ^ COnjUnCti ° n With ,he P-e-nnvention. 
with the present invention W ° ,0rcalcu,a,in 9 a " '"tensity profile and intensity gradient in accordance 

S RG n 0 *«„s™Ltf^ n ' "V"™ «««•«" ■» Present invention. 
DETAILED DESCRIPTION OF THE INVENTION 

K^^^^^^ of explanation, numerous specific details are se, forth in order 
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the intensity profile of the pupil is calculated using the kernel A[ijJ and the vectors corresponding to the mask parameters. 
In step S412, a simulated aerial image of the mask is produced. 

[0054] Prior to discussing the operation of a simulation system in accordance with the exemplary embodiment of the 
present invention in further detail, a brief but concise discussion of the mathematical underpinnings will be provided. 
5 Specifically, the spectral basis, including its analytical representation, its basic properties, and its application with re- 
gards to an indefinite integral will be discussed. Further, the Fourier transform ("FT') of the basis, including its analytical 
representation, and its application with regards to an indefinite integral will be discussed. Still further, miscellaneous 
provisions including orthogonality under a FT will be discussed. 

'0 SPECTRAL BASIS 

Analytical Representation 

[0055] A mathematical function provides a mapping of a continuous set of input values to a continuous set of output 
is values. This mapping may be generated in a number of ways. One way is the use of an analytical expression, such, 
as sin(x), or log(x). Another way is by tabulation and interpolation of intermediate values. Yet another way is to represent 
a function in terms of a linear superposition of a set of known basis functions. For an orthogonal set of basis functions 
such a representation is called a spectral representation. 

[0056] The spectral basis considered in the following discussion is represented as: 
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25 which may be factored into one-dimensional contributions. In equation (1), r = a (r-r 0 ) for sampling at a point r 0 , and 



0l [a] lr> = (a/^)(-^) Aa expfr / 2 ]a/a/exp[-? 2 ] (1) 

'x 'y 

dimensional contributions. In equation (1), 
X. a - a x + ay The one-dimensional basis along the x-dimension follows as: 

i 

0l [ax , l*> = («/VS) 2 {-h aX exp[x 2 / 2]d p exp [-x 2 ] (2) 

'x 

Basic Properties 
[0057] 

1. Recurrence 

^[a^") = ^rajl^a,^ Jlr> (3 ) 

2. Derivative 

^V >= I a ^[-i/j ,r >*< 1 1[ a+ i / ) ,r > (4) 

3. Normalization 

The basis is orthogonal with normalization constants: 



a.!.a.. f 



(?,(r)'7,(r))= with a=(a Xj a y ) 



x-"y 



(5) 



Indefinite Integral 

[0058] The indefinite integral over the spectral basis may be related to an incomplete Gamma function. For Manhat- 
tan-style layout designs, the integration boundaries are along the Cartesian coordinate axes and permit a factorization 
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^independent one-dimensional integrations. The following def.nition of the incomplete Gamma-function Y (a : x) 



Let 



JO 



•■= \:'dx(rj \ x ) 

• l"xl I 



[0059] Recurrence 



(6) 



(7) 



Let x = ctx and 
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[0060] Note that J (ajt m| (-i) =- (-1 )/» ♦ a (i) The above jn , egra| may thgn fae eva|ua{ed 

l [a,(>fo^l)=J [a ,o ) (il)-J [ a.o,(^o) 



(8) 



as 



(9) 
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J|>-l.»] = 
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[0061 ] The Fourier transform of the spectral basis may be expressed as: 



<11fall * = aiy) XaeXp[ ^ /21d ^^ expf ^ 2 J 



Indefinite Integral 



(11) 



(12) 



» [0062] Similar to equation (3), the indefinite integral over the FT of the spectral basis may also be related to an 

Z sZ Ti^ZT^ T t0 CircU,ara P erture °< < he ^ however the two-dimenJLa! in^egla fo c nnot 
dLrln V, . onc-d,mcns,onal products. Tho present invention uses a hybrid approach, whereby the ono- 

o'er hTc ^ aXiS ' f ° r 6XamP,e ,he X " aXiS ' iS P erformed ^ically s -i'ar to equaZ (3 On the 
s.ders the one-dimens.onal analytical mtegration of the FT of the spectral basis. Let 
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Recurrence 

15 



'w(*o.*i) ■= \ll d v{v { *A<l) (13) 

= <* Jl^v^ilq) (14) 

a 

= a i *- i 9 \ - a j *. 7 /fl/ ; 9 ) ( 15) 



[0063] Let k= k/aand 



■C,<*) 4^}{^J J" * e *P t 2 ^]a; exp [-r] (]6) 



then the above integral may be evaluated as 

25 



*[a)( k oA)= ^a<)](*l)-^a.O]<*o) (17) 

A recurrence for J* [a mj fly follows via 



30 



i ^(0 m exp[r 2 /2] 5f5 rexp[-^] 



(18) 



35 

as: 



J ^1™1 " <- i )Jf J [a,m + 1) + nvJ [a.m-1]-Pk m {il [a] lk>+8 raj0 p<ii M IO>J (19) 

[0064] This relation is identical to equation (10), with 0 = 1/« and an additional factor of (-/) per recurrence step. 
Initialization 

J5 [°° 651 The ab ° Ve seed - e,eme "ts J»io M are similar to the previous section and may be evaluated similar to equation 

l °' m) "l^J Y( -2-"'2 )=a Vm| (20) 

[0066] With the above refreshing review, the operation of a simulation system in accordance with the exemplary 
55 Inn Mh the present invention will now be provided in further detail. Table 1 , below, provides notations for symbote 

and their respective descnpt,ons that are used throughout the description of the present invention 
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Symbol 

r= (Ky) 
k = (k Xi k) 

, m,i 
<p 



M 

A 

s 



Description 

2-D position vector fixed z-plane " ~ ~ ~~ "~ 
2-D wave vector fixed z -plane 
plane wave entering mask 

scalar field at mask exit plane at normal incidence. 0 = 0 
scalar field at pupil entrance plane, general incidence angle 
scalar field at pupil exit plane, general incidence angle 
scalar field at wafer entrance plane, general incidence angle 
Mask transfer matrix 
Pupil transfer matrix 
Magnification factor (0.2, 0.25, etc) 
Scaled coordinate system at wafer plane 
Scaled coordinate system at pupil plane 
Partial Coherence 



(ol/) = i,(r) 



(21) 
(22) 
(23) 



35 



40 



<Mr> = </1/c> 
(rik) =(ri-k) 
(rlk-k') = <k'lr><rlk) 

[0068] Completeness of an orthogonal normalized basis £ implies 

J i = 1 

[0069] Spectral expansions used in the fo.lowing satisfy the above comp.eteness relation. 
Lens Model 



55 



(24) 
(25) 
(26) 

(27) 



Kb - T 



(28) 



[0071] Further, the apertures may be defined as (w= substrate, s = Source): 
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NA W = n, sin a w = sin «„, (29 ) 

NA S = n, sin a s = sin a s (30 ) 
_ NA S 

° = NA„ (31) 
[0072] The maximum wave-vector transmitted by a pupil of aperture NA thus follows as 

k NA =d, sinu„= d,NA (32) 
[0073] A magnification factor M may be defined as the ratio of the back- vs. the front-focal length of the objective lens: 

of, 

M = d 0 (33) 

SLt W t h , ere d °,T d , d1 <; 0rrespond 10 the focal ,en 9 tns identified in Fig. 1 A. The substrate plane tf, may be adjusted 
such that it lies in the focal plane of the lens : thereby providing: ujufiea 

1 _ 1 1 

7 = do + 5; (34) 

*» EX, ^ ,hem ^ tical model is ^pressed in terms of coordinates k s , r, k, V The numerical model however 

utthze. re-scaled coordinates indicated by "hat" notation. The scaled coordinates relate to the normalized coordinates 
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r m := - Mr m Mask plane 



C -—If. en A U ^ 



(35) 



s p '■= k 0 sin t> = — Sp Pupil plane (35) 



(37) 



a * s _ k 0 sin 9 S 



^ — Source plane (38) 



55 



The current lens model in accordance win the present invention assumes d s = d a 
Djffractive Elements and Lenses 

E, J* 6 . P r °P aaatio " ° f from the exit plane of a diffraction element, such as for example a mask, (+> Ir) to a 

H^ZT 7? P ' k) ^ diS,anCe maV be 6XpreSSed in lerms of lhe frea - s P a - Greens funct on of The 

Helmhottz equation ulGI_0(r= (x,y,z) in three dimensions) as 
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(r\<S>'') =ik 0 fdf(hG\Oir>W x )=ik 0 (r\G\V x ) (39) 

with 



-'k 0 lrVI 



4n\r-H (40) 



Mfr'^ko-K-^r-O 2 (41) 
[0078, Assume a far-field approximation .r - , , . d , and approximate by keeping to first order of a Tayior expansion 




(42) 



<rlG lr'> = " 



4wd ( . (43) 
Lens Element 

[0080] The propagation of light from the input plane to the exit Dlan* nf a i MC 

of the lens propagator L as P S element ma V De expressed in terms 

<rl<D''> = Jcfr'<rlUr'><r'l«), x >=<rim x ) (44) 
where the lens-propagator in the paraxial approximation may be expressed as: 

(nut) = e*' 5(r-nr(r) (45) 
P-^ Present invention may be impiemented using any 

operable to input lithographic projection aoDaratus D *n>m«*L m~ . . P y P erformed b V an y device 

inspection tools, and ^rnalmemo^ ™* 
operable to input the mask parameters Similar to steo sua 1 i.w, P ? V b * perfomied b V an V device 
inspection tools, and external memories comSnhntS? n ° n " ,,m,t,n9 e * a ™P>es include stepper/scanners, mask 
both thelithographicpr^ a -"ft* device may be used to input 

^1:™ J - step SS0, an area of the pupil of the 

gr,a - Al sle P S504 ' lhe illumination profile of each discrete point 
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w 



20 



25 



35 



45 



55 



corresponding to the sample area is determined, and given a numerical weight according to the relative intensity of 
the source at the sampled location. At step S506. an array of seed-elements for the semi-indefinite integrals calculated 
m equation (16) of each d ISC rete point within a pupil grid corresponding to the sample area is tabulated Each seed- 
a? ' S a " ' ncomplete 9 amma function J [o.m). such as is described above with respect to equation (1 9) and equation 
(do). At step S508. the pupil model is determined for the simulation method. If the pupil model is for a zero defocus 
pupiL i.e., for a lens having no aberration, is perfectly in focus, and/or has no phase distortion in the pupil plane then 
step S510 zero may be performed. If the pupil model has optical path differences, has aberrations is not perfectly in 
focus, and/or has phase distortion in the pupil plane, then step S510 optl may be performed. Further, if the pupil model 
accounts for diffusion of the photoactive compound in the resist after image exposure, then step S510 rfm may be 
perfonmed. In each one of steps S510 zero , S510 opd , and S510 dm , an array in 7, shifted by s, of p a ( S/ ) is calculated, 
wherein each entry is the value of orthogonal projection of the pupil onto the basis set indexed by 'a 1 At step S512 
the weights of each respective discrete point in the illuminator profile from step S504 are combined with the corre- 
sponding orthogonal projection on the basis set as calculated in any one of steps S510 2ero , S510 Dnd and S510h« 
The solution to the combination performed at step S512 renders the kernel A[ij] at step S514 

[0083] The imaging of a simulation system in accordance with the exemplary embodiment of the present invention 
will now be descnbed with respect to coherent imaging, incoherent imaging with the Hopkins model, and with respect 
to the basis set decomposition. 



COHERENT IMAGING 
Pupil Plane 



[0084] The field at the mask exit plane ♦"•'to given as the modulation of the incident field $ mJ by the mask transfer 
matrix T For an ,nf initcsimally th in mask lay#, Tis effectively diagonal (For a basis that satisfied I -i,> = a diagonal 
matrix <rl 71 r) implies that ft I r/n'>=/ft-iy)): J 



<^j x lr> = Jdr- lr->(r' I T lr>5(r, r") = 1 71r> (46) 
30 The field at the pupil entrance plane follows via equation (39) as 



<*fi]V = 'VCf IG ° ls P> (47) 



= K 0 <*w I TG° ls p > (48) 



Substrate Plane 



[0085] The propagated image at the substrate plane follows according to equation (44) as the field propaqated 
through vacuum segment 1 from the pupil-plane exit field <b px - Hip 

(si 



<4»r s ]'^>='V*fsriG 1 lr w > (49) 

H, J T h PUPil ; P 'f n f 6 f fie ' d itSe " iS thS PUP " plane entrance field O'isl Propagated by the lens transfer matrix 
{k I U I X). The substrate plane field thus follows as 



so f*":'ir \ = it./*?' n 



<* ls i'U=ik 0 <^ sl l IL 1 G 1 Ir,,) (50) 

T T ' enS tranSf6r matriX iS assumed to be Agonal in Fourierspace, (ML 1 !*■)= <ML' l*)6(x./r). Using equation 
(47) the substrate plane field may be expressed as a function of the field at the mask input plane i 



! as 



«fr s i i|r w=-ko<Ci i|TGOi - lGi|r v») (5D 
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Source Plane 
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incidence as: 



to 



15 



(52) 



sucH p la U n n e e ; a ve ilfa S^SS f 3 " ane — - norma, incidence. The imaginary part of 

0 reflecting the J^SSSSSi SSSS KST ^ ^ ^ iS ^ ^ ~ * 



aperture: 



[0] 



20 



25 



30 



35 



(53) 



[0090] Using the above illumination, the substrate plane field f 



' for P°' nt source located at /c s then follows as: 

<^w>=-w(O s )k* /dr m e iksrm (r m ITG°L 1 G 1 lr w > 
FfHunhofer approximation 

,2 



(54) 



(<l> t «jlr w )=-"<»,)k 0 Jdr^e' ^(r m l 71^<rJ G °, Sp >< Sp | L i, Sp><Sp , G ' lf , 

= "^?d^; JdrmdSp e ^ k S r -^< r n. ,Tlr m>PYs p > 



(55) 



(56) 



40 



x exp 



-ik, 



2d 0 2f, 2d. 



[0092] The argument of the exponential may be evaluated vi 



via: 



(57) 



45 



50 



55 



(58) 



^t^^uS^ zt'zn^s^ fa? fan out when squarina the fie,d to ° btain the ^ 

the image at the mask plane r ~ Wr the sub ? r ™ * " the SUbStra,e plane a PP™™tely , 

P ' * wr m. the sub-term fr 0 can be eliminated for the same argument. 



' mirrors 



Defocus Term 



[0094] The center term in the above described exponential relates 
focused system, ihis lerrr " 
a Taylor expansion yields 



recused sy Sl em, ihis term disappears V<^rt7nVto eauaiten 7™*e° f0cuscondition system. Fora perfectly 
a Tavinr »,n a „c inn HP according to equat.on (34). For a system under small defocus of magnitude £ 
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1 1 



d,+ ; d, d, 2 



(59) 



wherein a phase term accounting for defocus thus follows as: 



10 



4" = k -i-c 
°2d; 



(60) 



[0095] A more accurate defocus model follows by calculating the true optical path difference ("OPD") of the Gaussian 
'5 reference sphere: 



20 



^(c)=koCd-cr x ) 



(61) 
(62) 



25 Rest of Terms 



[0096] The first term in equation (58) contains the main imaging information and corresponds to a Fourier transfor- 
mation of the system. Introducing the scaled coordinates equation (35) this term follows as: 



30 



35 



X = ik s 

0 p 



d d 

\ \ oj 



(63) 



40 



is ( r -f ) 



(64) 
(65) 



45 



[0097] Defining 



so 



2d, J 



c = P(s p )exp 



2*„ 



(66) 



the lens system equation (55) then follows as 



55 
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(67) 

Define ^mOx^m^xSplPlS^^fls,), then 



ro 



(<;;i^)=w(^)(s $ |fp| rw ) (68) 

'5 INCOHERENT IMAGING: HOPKINS MODEL 

dependent point sources S : substrate plane follows as a weighted summation over all in- 



20 



25 



30 



'^"H'SM' (69) 

= K(^l°ni s ) W (s s )^ s ,7P|rj (70) 

- I * ! *'}(*' I TH,) I 7- 1 *X* \P\rJj dkdk> 



(71) 



Let 

35 



*o then 



(72) 



Kr w) = Jdqdq'dk, (rJP|q' + g s >t(q')w(g s ) 2 f(q)< q+ g| P | rJ 
since {K\P\K> is diagonal define P(K) :=(MPIk): 



(73) 
(74) 



50 



« r .)=H"q'.( q ') t .(,)Jd k ,(r.| q . + S>)P . (<I . + §> ( l<)!(q+5i|r _ )I>(q+gi) (?5) 
-"** W *'W('-.l' , -1>I'*.P-(»'«.W.) , Pfc«,) (76) 



15 



EP1 202 119 A1 

Define the Transmission Cross Coefficient (TCC) function T(q\q) 

T(q',q)= /dk s P*(q'4-S s )w(S s ) 2 P(q + g s ) (77) 

then 

'° 7 (0 = \ d * d <l\\A<!)tW)c{q\ q )t m ( q )( q \r w ) (78) 

Assume a factorization exists such that 

15 

T(q',q) = X^(q')A a i lj (q) (79) 
50 where ^q) forms an orthogonal basis. Combine equation (77) and equation (79): 

Z<(q>^i^ (80) 

55 0 

[0099] Let N, = (n /T1 ,) be Z.2-norm of r,,. Left- and right multiply equation (80) with the orthogonal basis n, , followed 
by an integration over the arguments: J 

Vfpr/ dq d q'*.1,(q , )P*(q , +6.)w (s s ) 2 P(q + S s )^(q) (81) 
[0100] An orthogonal pupil projection coefficient^^ may then be defined as: 

Pi(k) ^JdqPfq+C)!,:^) (82) 



30 



35 



40 



^fdtfpftritf . K) (83) 



[0101] in accordance with the present invention, p X *) may be pre-calculated for any possible pupil configuration i 
e.. cons.denng non-limit.ng parameters such as NA, focus, aberration, etc. The radial part of the P(k) is typically a 
c.rcular aperture. The phase part describes aberrations of the lens field. 

[0102] Using the orthogonal pupil projection, the coefficient matrix, or kernel. A ff of equation (82) follows as: 



A, = /dk w(k) Pj *(k) Pj (k) (84) 



Note that A,y is Hermitean. 



[0 03] Stop S506 may bo performed in accordance with abovc-describod equations (20) and (1 9). and as illustrated 
wrfh respect to FIG. 7. In a first stop S702, tho column a = 0 is initialized in accordance with equations (1 1 ) or (20) 
wherein each entry in the column is seeded with an incomplete gamma function. In the next step S704. the values of 

in° hT^oi = rMT 393 ^ t0 th ° C °' Umn l ° ,hS ri9h ' in accordance Rations (10) or (19), wherein each entry 
in the column is filled according to recurrence. After repeating step S704 until row m = 0 is full the resulting indefinite 
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5 Algorithms: Aerial Image 
a. Pupil 



[0105, Combining equation (82) and eouahon (66). the orthogonal pupi. projection coefficient fol.ows as 



1 r. r /i' 



^^)=^-/ d 4P'(q)e X p^c 



20 



Approximation 1: Aberration Model 

I ^ m ° dU,alin9 ^ PhaSe M ™ — P"P« opening. Normal the 



f:=i 

k 0 (86) 

25 



using = 2 * 1 the phase error K caused by defocus can be 



expressed as: 



• 4 (f. t ) :.JL ?= i t< ,,( 2 , +z J 



JO 

(87) 



35 



* ( '' !) 'I*. (88) 

40 

[0108] The orthogona. pupil projection coefficient then follows by numerical integration of 

P - (I ^> ^^«pNf4n;(q-s,) m 

[0109] For terms that have on.y radial dependence <o :=.$,) the integration may be simp,ified as: 

50 

P.(«..S) ^rqdqexpMf.^jf^n^q-sJ (90) 



55 Approximation 2: Zero Defocus 

ations, the pupil fui 

1 circular integration 
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P.(s,.s) = t?t \dqexp 



nl(q-s s ) (91) 



[01 1 1 ] For a lens system that is in-focus (C = 0) equation (85) reduces to a convolution of the pupil aperture with the 
spectral basis and p (q„s s ) 



to 



(92) 



>5 [0112J The circular aperture may be approximated by a superposition of rectangular domains B. as illustrated in FIG. 
3. FIG. 10, on the other hand, illustrates the improved method of the present invention, wherein for each domain B 
the pupil convolution evaluates to 



25 
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20 Pa(S s )=P a ,(S s<x ) Pay '(S s>y ) (93) 



(94) 



[01 13J An efficient method for evaluating the above one-dimensional integral, will now be described with respect to 
FIG. 11. 

[0114] From equation (17) the zero defocus case may be expressed as 

30 

Pa(S s )=^(<.<)^(^y.dq;) (95 ) 
[01 15] Sampling over a finite grid /c s , a finite approximation to A l7 follows as- 

35 v 

A.. = A^ZwCiJ^t'JpyfeJ (96) 



[0116] Note that r, , (A), and thus p.(s s ). are purely real for even indices and purely imaginary for odd indices The 
.mage amplitude is evaluated using only the terms 2Re Thus, for in-focus conditions all entries in A s that are 
purely .maginary may be disregarded, as illustrated in FIG. 1 1 A, as sub-block 'D'. The remaining diagonal A and sub- 
blocks B and C may be organized in software as illustrated in FIG. 1 1 B. 

Approximation 3: First Order Defocus 

[0117] A first order approximation towards defocus aberration follows by calculating a series expansion of Pa(s ) 
around C, = 0 up to the first order lean: 1 ** 



dc 



(97) 



[01 18] The first term is given in equation (95). The second term may be evaluated as: 
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^tt -s s ) (98) 



45 



50 



2k N i uqq ^.W-sj (99) 

0 t o 



[01 19] For C = 0 this term follows as: 

as P .(s,. ? = o) ='_i__j d q^ T):( ^_ s j 

[0120] A first and second order correction to the image kernel follow as 

K = A ' i .?-(«.)[P;(i.)9 f P j (*.) + (^P;(*.))P i (8 i )] (100) 

A " - (ioi) 

He^lan^ ST"^ ^ 6 ' ementS PUre,y •"«- while A«K is 

coherent vector o ^^^^^^L h t "t" ° f ** M * y ' e,dS * bilinear f °™ ^ • 

even-odd and odieven eZents 2 of e ^evance while for AV 0 £ IT'* * * 9 " ^ F °' ^ thUS *• 
to the image intensity. relevance while for 4« # only the even-even and odd-odd elements contribute 

*> [0122] The first order correction to the image intensity then takes the form: 



lM = _S( r '»>-(»-. |r > + ("»-> A r i <«-l^ 



(102) 



E ss^f 6 ^ ?»■«** calculate a number of image planes throughout the 

distortion in the outer pupi mg vs 2 S me H^orde °T ***** « Sma " (comparing the phase 

*> behavior through a t^I^SS'JSSS^ approx,mat,on ^ a efficiently accurate mode, of the defocus 

Algorithm 4: Improved OPD integration scheme 

[0124] For a general solution of equation (91 ). the following steps may be used: 

1 . Discrete the domain of the pupil over a finite grid with indices k x , k y s x , s y . 

2. Build a matrix of OPD values in the pupil plane ^ : . c/^expW*,, kyQ] 

3. Build vector of basis function values in first dimension: -4(a x , q x ) := ^ " ( q x) 

4. Calculate the following vector: 



55 
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C ( k »> a -)[ s ,] = E M P (k,,k y ,^)A(k y -s y ,a v ) 

(103) 

5. Use the above matrix of vectors to combine all calculations in the second dimension: 

l" n 

(104) 

Incoherent Imaging: Hopkins Model with Diffusion 

[0125] In the discussion immediately beiow : the Hopkins model is modified to include the effects of photoactive 
compound diffusion in the resist. An isotropic Gaussian diffusion model is assumed. Diffusion is modeled as convolution 
with a Gaussian kernel: 



2ko 
with 

a = JlDt 



2 e (105) 



(106) 



[0126] Assuming the photoactive compound concentration is linear in /(/). diffusion may be approximated by replacinq 
equation (73) with: . y K y 

<r m i-Hr m \Zl (107) 
[0127] Comparing the above to equation (70), this is equivalent to substitute 

P^ZP (108) 
[0128] The intensity in equation (74) then follows as: 

l(r w ) = / dqdq't(q')r(q) /dk s <r w |ZP|q ! + a s >w(a s ) 2 <q + g s |PZ| r w > (109) 
[0129] The bracketed term inside evaluates to: 



20 



\>\.\ZP\q' + 



(110) 
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V + i,} = J*A'(r.S*X*|Z|Ar')(*1/»|^ + ji ) 

= ldk(r^)( k \ Z \ q - + s s )p{g + s 5 ) (H1) 

[0130] If ( r/2l O Mr- then z is diagonal in k-space. Thus: 

<rjzp\ q ;-$ s) = a s)Z(£J , A s)p(<; , ( ^ ^ 



20 p i(k) = ^/dqZ(q + k) P ( ( ^ 

[0132] As relating to FIG. 6: 
25 Point Decomposition 
[0133] Define: 

30 < r w , 9i> = /dk<r w IkXtlkXkliij) 

[0134] Then equation (78) can be expressed as: 



(113) 



(114) 



35 



40 



45 



<ffl -I rj amounts to a convolution integral in real space: 

<9/l i) =ldr , x [j \ r-t-)(t\ry 
[0135] Decompose the mask transmission function into a disjoint set 



(115) 



(116) 



(117) 



50 



55 1 P = V/1p 



v (118) 



and thus 
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(119) 



Reduction of Field Integral 

[0137] The above expressions reduce the mask transmission function to a summation over an indicator function 1 
Furthermore, for Manhattan geometry's, this indicator function factors into the cartesian product 1 P = 1 Px • 1 p The 
next step is to simplify the integration in equation (117) utilizing the above reduction of the mask transmission function. 
[0138] Factorization of the indicator function is matched by the choice of the spectral representation as discussed 
above with respect to basis set decomposition. The basis r), itself also factors into a Cartesian product n.„ = ti * • ti y 
• This permits decomposition of the two-dimensional field integral of equation (117) into a product of two^ne-^imen* 
sional field integrals. Define: 



/„ (dp x °, dp x ] ) := Jdx ' (< | x - x')(l p J x ') 



(120) 



(121) 



- Ji 



\>?:, i*-; 



(122) 



[0139] Where 



dp x := x - p : 



x,0 



lower boundary 



(123) 



dp x := x-p x A upper boundary 



(124) 



[0140] Using the Gamma-function factorization as discussed above with respect to the spectral basis, the above 
integration may be reduced to a recurrence over the semi-indefinite integrals J a (x) via equations (4), (9) and (10). 
[0141] The above expression equation (117) then simplifies to: 



(125) 



[0142] A particular advantage of this representation follows when several layout features line up along a Cartesian 
dimension. Assume for example q e n, to be a subset of features that are vertically stacked up and have identical 
transm.ss.on t, M . Then for all q the values of q x0 = cq and q x , = q are constant. The above expression then reduces to 



C^'lr) = ^l ni (dc:,dc , K )Xl r Jdq°,dq;.) 

qeflj 



(126) 



intensity Calibration 

[0143] For a fully transparent mask : I p 1 p (r)= 1, f p =i and equation (11 7) yields: 
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Vn v^ ; ' (128) 

Dark Field Reticle (binary) 

* P y (129) 

75 Clear Field Reticle (binary) 



20 



25 [0146] Thus: 



Clear Field Mask (binary attenuated phase shift) 



35 



45 



(132) 

simulation domain: 



40 

[0148] Note again that l+ p and 1 - p completely cover the 



(133) 



55 



[0149] Thus 



(134) 
(135) 
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Clear Field Reticle (ternary attenuated phase shift) 



[01 50] A ternary attPSM clear field reticle consists of three transmission domains: a 0° phase fully transmissive clear 
field region 1 a 180° phase shifted region 1 with transmission a. and a chrome blocking region 0 o . A full partition of 
5 the mask area follows by the disjoint union of ail three domains: 



w 



p p' n 



(136) 



75 



[0151] The mask transmission function can be expressed as 

<?\r) - ZCto-tfX.i;. (137) 
p p 

= ] -I.K-^- p -T% (138) 

P P' q 

= Wa+DlX-X.O, (139) 

P q 

[01 52] The coherent field vector follows as 

*> f> = {g; n>-(a + 1)<of" lr- <pj*~~ /r> (140 , 

Dark Field Reticle (ternary attenuated phase shift) 

« [ ? 1 *? ,7 he d ° main partition of a temar V attPSM darkfield mask is identical to the domain partition of a ternary attPSM 
Clearfield mask. The mask transmission function however is expressed in terms of 1- and 1- : 



25 



40 



('lO = Li;0-)-*Zi;- (i4i) 

p p 
and thus the coherent field vectors follows as: 

tof**""*!')-^ !/>-•<«; i/> 



(142) 



[01 1 54] For an attPSM contact mask, the number of rectangles can be reduced through a different linear combination 
of 1* and V. For a contact at index p define an 'outer frame region' 1/that includes both the 0° and the 180° region as: 



1 p = 1 p +1 p (143) 



the mask transmission function then follows as 

55 
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<'"■> - ZW-al- (144) 

- IW-^-i;) (1 45) 
= Z(i+*)i;0-)-< (146) 

the coherent field vector can then be expressed as 

(tf^^UaXg^.^ (147) 
Fast implement ation of Zero and OPD Model 

^e^^ -i^ons in the 2 ero- 

m (x.y) space. As illustrated in FIG. 10 th feel o ^"^ to ^ fcu h ^«? ,n ^ 0 «^.ci ro utarp U p l Ww 
algorithm as illustrated in FIG. 1 0 is as follows Rrstt » T " Y ^ ^ me,h ° d f ° r e erf °™"° the 

intersecting a ray emanating from ^ ccntcNs M ? 7 T*™* °' 45 °' WhSrein the P ° int at the cire,e ™ 

across the x and y axis throug the center S he SL Th " °. th ° thr ° C P ° intS aCrOSS ,hc contor 

<rom the center at the current ang e is 2 ^22 AaafnT PWn, . at CirC ' e " B in,ersecti "9 » «f ^anating 
rectangles is determined Then the contribuZ of th 9 ^ymmetncal point is determined and the new set of 

is determined. contribution of the new rectangles to p[a] according to the equations (95) and (1 7) 

S^tSCSiSSJSlSr 3 a Thirbec h a e PUPi, ,h aS dSSCribed ^ re,erSnCe ,0 F ' G - 10 * ™<* -re 
is no. calculated, thus a large Amount of ca^II T " m6th ° d ° f RG - 1 °' eVery discrete P° rtion <* the pupil 
number of rectangles to cSClX^ SPeCifiCally ' " * 1 ° ^ 3 

SseorseTa:' 

for performing this step are known in the art Steo 2 hll !f eC,f,ed f ° rthG aeome tric sampling. Methods 

<° [01 59J As opposed to arranoinc nhJLT ',1 P ' however ' 13 Performed contrary to that of the prior art 

•he prior art systems the^Znt Z^ ^ZT^ 'V™** ^ °" the ^ 0f the ™ sk - " 

uses a one-dimensiona, «^7^r22^TSEJTliTr^r n ~^ 0ther Mb ' ,he pri ° r art method 

resulting in an array of tmJk^StS^ Sn ^^l TT™* ana ^ onhe geometrical indices, thereby 
dimens 10 na, array for the ifS^Sl?^.^!^ T m6,h ° d °' ^ pr6Sent invention uses a one- 
's in an array of one-dimensiona. tables one-d.mens.ona. array for the geometrical indices, thereby resu.ting 

scribed above wi.h respect to eouaZ mf «h Seed ' e ' ement ,S an '"complete gamma function J[a], such as is de- 
field exposure f^^JS^^^SS^' ^ ""J iniUa,i2ati ° n SleP S6 ° 8 Ca ' CU ' aleS lhe Clear " 
» tion (9). a ' Ue forsub sequent exposure dose calibration with equation (128) and equa- 

^LlToZZllTe T^j^T ' n r er, ° r ^ aerial »™""'" "«* pattern to besimu.ated, 
RG. 8, item 800 is a predefine proxTmi* ^wtZ " n0,VOVerlappina Wangles. As further illustrated in 
Plurality of samp.e features 804^2 %™ ?SZ ZT. ™ ^ °' ltCm 8 ° 2 ' 3 Samplc mask comprising a 

■ -e^ 

variatL o^^^^l^Z ?"™ T * """"" d ° mai " ^ the edge 
bas.s (g 0 . ) remarns .arger than or equal to <= aec . The boundary of this domain follows as: 
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2 2 12 



( f i.p«,x)-= =e (148) 



Thus 



r L p TO x= (149) 
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[0163] In essence, the proximity window defines the set of all rectangles that affect the sampled area. Specifically, 
in this example, as illustrated in FIG. 8, features 81 4-818 do not have an affect on the aerial image of features 804-808! 
Therefore, features 814-818 are not within the proximity window. 

[0164] Step S612 will now be described in more detail. Once the features within the sampling window have been 
decomposed into non-overlapping rectangles, the corresponding orthogonal projections of the features of the mask 
onto the basis set are calculated in accordance with equation (125) and equation (9). 

[0165] Step S614 will now be described in more detail. Once the orthogonal projections have been calculated in step 
S61 2, corrections may need to be applied, for example, depending on the mask. If the mask is a clear field mask, then 
equation (131) may be used to correct the orthogonal projections. Further, if the mask is an attenuated phase-shift 
mask, then equation ( 1 40) may be used to correct the orthogonal projections. Finally, if the mask is a ternary attenuated 
phase-shift mask, then equation (142) may be used to correct the orthogonal projections. Once the corrections have 
been applied in step S614, the final corrected orthogonal projection of the mask transmission function g[/] is attained. 

Intensity Summation 

[0166] Returning back to FIG. 4, once steps S404 and S408 have been performed, the image intensity distribution 
may be calculated at step S410. Step S410 will now be further explained with respect to FIG. 9. As illustrated in FIG 
9, at step S906, the A[ij] from step S902 is multiplied with the g[/] from step S904, in a bilinear form. This calculation 
provides the image intensity l(r) at point r. Furthermore, the intensity gradient at point rmay be calculated at step S908. 
As illustrated in step S912, the intensity gradient may be measured with respect to the gradient in either the x or y 
coordinate axis. 

[01 67] Specifically, the image intensity l(r) at a particular sampling point, such as in step S906 . can then be calculated 
as 



'00 = Z(r\9,)A nm (g m \r) ( i 50 ) 



Mr) = 2Re£^<r|ff.)(f(^ 1 |r)-(ff : ^ IiJ |r)) (151) 



[0168] A numerical implementation may take advantage of the hermiticity of A nm . The basis r\£i) is real. For binary 
masks, as well as 180° phase shift masks, t p is additionally real. Thus, the number of summation terms may be cut 
approximately in half; 



= HM 



rig 



•/V 



» 



(153) 



Dose Calibration 

55 

[01 69] For a computational implementation a calibration may be used against a clear field area for absolute exposure 
calculations. For example, using equation (128), the intensity of the clear-field exposed area then follows according to 
equation (150) as: 
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L + W\S't) (153) 
[0170] Therefore, given an exposure dose D, the local intensity sampled at point rear, then be calculated as: 



D ( r ) = ^Z^Ms,)<S j \r) (154 ) 
Th^inr^n 0 ^ 03 '? 131 ^ ' WhiCh 13 ,hS imenSity calculated b V the non-normalized algorithm for 100% exposure 

]J JESZ^X te the systems 10 9,ve an exposure dose va,ue ° as specified by the user - a 

!°ImL c T i h n e , S H S,em a " d meth ° d ° f thS preSent inV6nti0n ma * be incorporated into many applications. Non-limitinq 

SSi^T T 6n f e " ^ 01 an M ° del ° PC Ca ' CUlati0n; Simula,i0 " e "9 jne for P-dicting an" 
Sna Lnd ! Peelers, ««ch as NA, illumination, etc., for drawn layout features; simu.ation engine for pre- 

SZ T a 2 T I"' f6atUreS ' SUCh 98 S6rifS ' SCatterin9 barS ' etC ' 10 enhance the PnntabiMty and 
process ^ I sfmX on T h" en9 ' ne * ^ P"""" 9 ° f defeCtS ° n 9 maSk for a s P ecified ideation 

optimize ,he ba,ance - - a — -p— * - «- 

mlnuL,nI h °f 9h t SPCC . ifi ! rC,CrCnC ° b ° mad ° th ° tCXt 10 th0 uso of 'i«h°9raphic projection apparatus in the 
manufacture of integrated circuits, it should be explicitly understood that such apparatus has many other possible 

oattems f ^P' 6, il may be em Ployed in the manufacture of integrated optical systems, gu Wance Ind detec u^n 
patterns for magnet.c doma.n memories liquid-crystal display panels, thin-film magnetic heads etc The skilled artisan 
Z1Z T l l n 'I" C ° meXt ° f SUCh alternatK,e a PP |ica ti°ns. any use of the'terms "reticle" or «SS TtSZ 
f017^ VZT ^"9 -placed by the more general terms "mask" or "substrate", respective.y. 
used to nlnZ 2 ™h 3 k tSXt h3S concentrated ° n lithographic apparatus and methods whereby a mask is 
T n ?' °k lat '° n b6am 6nterin9 the Potion system, it should be noted that the invention presented here 

toTJn r n br ° ader C ° nteXt ° f '^graphic apparatus and methods employing generic kerning mean? 
be u Sed to radiatl ° nbeam The te ™ "Pining means" as here employed refers broad.y t mean thTcan 

la^ nTZ a " : nC ° m : ng radiati ° n bGam With 3 patterned ^oss-section, corresponding to a pattern that is to be 
nLntn if 9 °' SUbStrat6: the t6rm " li9ht Valve " has also bee " used I" this context GeneraNy the said 

pattern will correspond to a particular functional layer in a device being created in the target portion such as an 39 
mtegrated cram or other device. Besides a mask (whether transmissive or reflective) on a mas ^tab^e suchoatternin! 
means include the following exemplary embodiments: ' P attermn 9 

L 1 l^^° 9 , ra M mable T" array - ° ne eXample 0f Such a device is a matrix-addressable surface having a viscoe- 
last.c control .ayer and a reflective surface. The basic principle behind such an apparatus is that (for exampte) 

fncSnnlT 38 1 ref,eCtiVe SUrfaCe refl6Ct inddent N9ht 35 diffracted wha reas -addressed areasS 
aCt6d li9m USin9 an a PP r °P riate m ^ the said undirected light can be filtered out of the 
[TZZT 0 1T 9 °T dWraCted li9ht b6hind; in thiS manner ' ,he beam becomes P a "emed accordi g to 
Irrav 9 T °' maiMe ^ surface. An alternative embodiment of a programmable mi?ror 

TIuLhT, y , H a 7 arra ? gementoniny ^ 

able such SSr" h " by , emP ' 0ying Pi -° eleclric ^ans. O-e again, the mirrors are matrix'addS 
mi ™ th addreSSed mirr0rs w '" ref,ect an ^ming radiation beam in a different direction to unaddressed 

mL on qU ' r matr L X addreSSi "9 can be P erformed ^ing suitable electronic/computer means. More tnfor 

and US S 53 T o r y H :, h r rSferred l ° Can 66 9leaned ' '° r 6Xamp,e ' from Uni,ed P «ents US 5.296 89 

b^ronce 3 Pat ° nt applications WO 98/38597 and WO 98/33096, which are incorporated heroin 

wni^ P il° 9rammab ! e ^ D afray An eXamp ' e ° f SUCh 3 construction is given in United States Patent US 5 229 872 
which is incorporated herein by reference. a,«»,o/*, 

[0174] Any reference to the term "mask" in the Claims and/or Description should be interpreted as encompassing 
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the term "patterning means" as hereabove described. 

[0175] In summary, although certain specific embodiments of the present invention have been disclosed, the present 
invention may be embodied in other forms without departing from the spirit or essential characteristics thereof The 
present embodiments are therefor to be considered in all respects as illustrative and not restrictive, the scope of the 
invention being indicated by the appended claims, and all changes that come within the meaning and range of equiv- 
alency of the claims are therefore intended to be embraced therein. 



Claims 



2. 



A method for simulating an aerial image projected from an optical system, said optical system including a pupil 
and a mask plane, the method comprising the steps of: 

providing a mask to the mask plane; 
obtaining parameters for said optical system; 

calculating a kernel based on an orthogonal pupil projection of said parameters of said optical system onto a 
basis set; 

obtaining parameters of said mask; 

calculating a vector based on an orthogonal mask projection of said parameters of said mask onto a basis set- 
calculating a field intensity distribution using said kernel and said vector and 
obtaining aerial image data from said field intensity distribution. 

The method of claim 1 , wherein said parameters for said optical system include aberrations. 

3. The method of claim 1 or 2, wherein said step of calculating a kernel corresponding to said parameters of said 
optical system includes the step of tabulating an array of incomplete gamma functions corresponding to respective 
points in said pupil of said optical system. 

4. The method of claim 1 or 2, wherein said step of calculating a kernel corresponding to said parameters of said 
optical system further includes the step of tabulating an array of orthogonal pupil projection coefficients corre- 
sponding to respective points in said pupil of said optical system, 

wherein said optical system is in-focus. 

5. The method of claim 1 or 2, wherein said step of calculating a kernel corresponding to said parameters of said 
optical system further includes the step of tabulating an array of orthogonal pupil projection coefficients corre- 
sponding to respective points in said pupil of said optical system, 

wherein said optical system is either not in-focus or has aberrations. 

6. The method of claim 1 or 2, wherein said step of calculating a kernel corresponding to said parameters of said 
optical system further includes the step of tabulating an array of orthogonal pupil projection coefficients corre- 
sponding to respective points in said pupil of said optical system, 

wherein said optical system accounts for effects of photoactive compound diffusion in a resist layer present 
on a substrate onto which an image of the mask is to be projected. 

7. The method of claim 4, 5 or 6, further comprising the step of combining sample weights of an illuminator profile of 
said optical system with said array of orthogonal pupil projection coefficients. 

8. The method of claim 1 , wherein said step of calculating a vector corresponding to said parameters of said mask 
further includes the step of specifying a proximity window within said mask for geometric sampling. 

9. The method of claim 1 , wherein said step of calculating a vector corresponding to said parameters of said mask 
further includes the step of decomposing a geometric pattern of said mask into a disjoint set of rectangles and 
tabulating an array of projections of the rectangles within a proximity window. 

10. The method of claim 9, wherein said step of calculating a vector corresponding to said parameters of said mask 
further includes the step of correcting the array of projections of the rectangles based on the type of mask. 

11. A method for simulating an aerial image projected from an optical system, said optical system including a pupil 
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and a mask plane, the method comprising the steps of: 

providing a mask to the mask plane: 
obtaining parameters for said optical system: 
obtaining parameters of said mask- 

0 ^Z^lS.C n ° m ^ P °' yn0mialS inC ' Uding P ° ,yn0mialS aSS ° Cia,ed ^ ^ ,ers 

^^X^SZ^- 00 ^ 00 funct,on associa,ed with said optical ba - d - - 

12 «ss srrpssassr an optical system ' said opticai sys,em 

a first parameter obtainer for obtaining parameters of said optical system- 

o^S"" SET keme ' b3Sed ° n - «"*" - - Pa-eters of said 

a !™ P ara ™ er , obtainer for obtaining parameters of a mask provided to said mask plane- 

:^;sz^T ,atin9 a vector based on an orth ° 9onai mask — "P— — - 

a third calculator for calculating a field intensity distribution using said kernel and said vector and 
an aenal ,mage obtainer for obtaining aerial image data from said field intensity distribution 

13. The processing device of claim 1 2, wherein the first, second, and third calculators are the same calculator. 

14 r„=nrp= ^^^a , =ssr an op,ical sys,em ' said opticai sys,em 

a mask within said mask plane; 

a first parameter obtainer for obtaining parameters of said optical system- 
a second parameter obtainer for obtaining parameters of said mask- 

w i f t l h rS ^ a '^ lat ° r . ,0r °T ° 90na " y projec,in 9 Polynomials, said polynomials including polynomials associated 
with said parameters of said optical system and polynomials associated with parametei 2 

15. A device manufacturing method comprising the steps of: 

llTti~ 

ZP^J^&gS ^ 10 " *»" Said ^' ^tem is simuiated 

cato«oTkTrlT h f ° r T ° PtiCal SySt6m ' Wh ' Ch C ° mpriSeS a P u P i! and a mask P'ane: 
omoSs set S " " 0rth09 ° nal PUPH Pr0jeCUOn ° f Paramet6rS * ■" ■*•»■ system 

obtaining parameters of said mask, provided at said mask plane 

calculate a vector based on an orthogonal mask projection of said parameters of said mask onto a basis 

calculating a field intensity distribution using said kernel and said vector and 
obtaining aerial image data from said field intensity distribution. 

55 16. A method of creating a data set representing a mask image, the method comprising the steps of: 
providing a data set representing a proposed mask image- 

simulating the aenal image of said proposed mask image using the method of any one of claims 1 to 1 1 or the 
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processing device of any one of claims 12 to 14; and 

adjusting or verifying said data set using the results of said step of simulating. 

A method of manufacturing a mask comprising manufacturing a mask using a data set created by the method of 
claim 16. 

A computer program comprising code means for causing a computer to operate to perform the steps of the method 
of any one of claims 1 to 1 1 . 
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Fig.lA. 
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Fig.2. 



Input 
parameters 


"^-202 


i 


f 




Approximate 
TCC 


204 


i 


r 




FT 

Reticle 


"~^206 




r 




^TCC 
FT 


^-208 



Fig.3. 















































































\ 




























\ 
































r- 




























-\ 






















































































































































rj 




























































1 


\ 






















) 


7 

















































32 



EP1 202 119 A1 



CM 



■■■■■ ■ 
■■■■■ ■ 
■■■■■ ■ 



■ 



D) 
CO 



.2? 

'l. 

< 



2 >> 

ft 

cu 



o 

5" 



o 



o 



2 

CO 

O 



CO 

O 3 

Ic "c5 

a S ™ 

£ c | 

= o J5 

LU 2 



3 — 
CO 

o 



00 

o 



^ CO 
(O Jn 
CO <B 

!i 

c co 
LU q_ 



CD 

o 



33 



EP 1 202 119 A1 



Fig. 5. 
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Fig.6. 
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Fig.9. 
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Fig. 11 A. 
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